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Abstract 

The large hierarchy between the neutrino mass scale and that of the other fermions seems 
to be unnatural from a theoretical point of view. Various strategies have been devised in order 
to generate naturally small values of neutrino masses. One of these techniques is neutrino mass 
generation at the loop level which requires a mechanism, e.g., a symmetry, to forbid the lower 
order contributions. Here, we study in detail the conditions on this type of symmetries. We 
put special emphasis on the discrete Z n symmetries as a simple example but our results can be 
also extended to more general groups. We find that regardless of the details of the symmetry, in 
certain cases the existence of a lower order contribution to neutrino masses can be determined 
by the topology of the diagrams with a given number of loops. Finally, we discuss the lepton 
flavor violating rare decays as well as (g — 2) M in this class of models, which generically appear 
at the one loop level. 
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1 Introduction 



Explaining nonzero but tiny neutrino masses is one of the most compelling open questions in 
modern physics. Various beyond Standard Model (SM) theories have been developed to address 
this question. The most famous mechanism for explaining neutrino masses is the standard (type I) 
seesaw mechanism [ I ]. In the standard and simplest realization of this mechanism, the smallness of 
neutrino masses is connected to the very large mass scale of new SM singlets (right-handed (RH) 
neutrinos). These new particles, being too heavy, cannot be produced at the LHC or any other 
man-made or natural environment, (maybe) except for the early universe, making a direct test 
of these models impossible. With the start of the LHC data release, it is more exciting to move 
towards models whose new particles are within the reach at the LHC. The smallness of the neutrino 
masses is not related anymore uniquely to the very heavy mass scale of the RH neutrinos, but 
requires additional suppressions, e.g., small Yukawa couplings, quasi-conserved lepton symmetries. 
A very interesting possibility is to forbid neutrino masses at tree level and have them generated 
at loop-level. The first proposals of this type of radiative neutrino mass models are the Zee model 
at one loop [2] and the Zee-Babu model at two loop [3]. 

All seesaw models lead to the effective dimension 5 operator, (HL)(HL). However in general, 
neutrino masses are not necessarily explained by this dimension 5 operator. Various other AL = 2 
operators can also give rise to neutrino mass. These operators have been classified in [I]. A 
particular class of these operators are H) m (H L)(H L) with more than one pair of Higgs fields 
attached to the corresponding diagram [5]. Recently, there has been a complete classification of one 
loop diagrams leading to the effective dimension 5 operator in [ ]. Obviously, the loop generated 
neutrino masses receive further corrections from renormalization group running, which have been 
studied for Ma's scotogenic model in [7]. 

Let us suppose that, thanks to a specific structure of the model, up to the n th loop-level, there 
is no contribution to the neutrino mass matrix. Of course, increasing the loop order will further 
suppress the neutrino mass. At the three-loop order, with M^ew ~ 100 GeV and couplings of 
order of 0.1, by dimensional analysis the neutrino mass will be in the range O ((g 2 /16t: 2 ) s Mnew) ~ 
O ((g 2 log(A/M NEW ) 2 /16ir 2 fM NEW ) ~ (0.01-1) eV with a cutoff scale A ~ (1-10) TeV. Such a 
high value of coupling and low mass scale is very interesting from a phenomenological point of view 
as it can lead to observable effects in colliders and indirect searches of new physics. Furthermore, 
the couplings leading to the neutrino mass can also in principle induce Lepton Flavor Violating 
(LFV) rare decays and a contribution to the anomalous magnetic moment of the muon. The 
construction of these radiative seesaw models often requires the introduction of an additional 
symmetry, G u , forbidding the tree- level contribution as well as contributions from lower loop 
orders. An interesting consequence of these symmetries is that they stabilize some of the new 
degrees of freedom [8] and these models can provide a suitable Dark Matter candidate [9-14]. 

In this paper, we will restrict ourselves to models which lead to the effective dimension 5 
operator. Hence, the left-handed lepton doublets L are the only fermions coupling to the new 
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particles, i.e. leptons act as a portal to the hidden sector. Employing a G u symmetry, we restrict 
the couplings to the form Cy = LSiFj, where Si (Fj) are new scalars (fermions) and forbid 
couplings of the form LHFj as well as LLSi. In this context, we will be general and will not 
restrict ourselves to the content of a specific model. Our aim is to outline general restrictions and 
no-go-theorems as a guide to build radiative neutrino mass models. We demand that all the SM 
particles are invariant under the new symmetry, G v , but some or all of the new particles transform 
under G u . In particular, we assume all new particles that couple directly to L as well as all the new 
neutral fermions carry a G u charge. Moreover, we assume that none of the new scalars receives a 
vacuum expectation value: that is, G u remains unbroken. As a result, the G v symmetry forbids 
a Dirac mass term for the SM neutrinos both at the tree level and at all orders of perturbation 
theory. The neutrino mass term should be therefore of Majorana type which in turn requires 
lepton number violation. 

For concreteness, we first consider Abelian Z n symmetries (i.e., we take G v = Z n ). We classify 
the emerging topologies up to three loop order and discuss the conditions on the Z n symmetry 
which forbid all lower loop orders. We then show that some of the results we find also hold valid 
for a !7(1) symmetry and more general symmetries G v in sec. 4. 

The paper is organized as follows. In sec. 2, we outline the general setting of the models 
discussed in the present paper and some general implications for neutrino masses. In sec. 3, we 
discuss the loop contributions to the neutrino mass matrix and show how the Z n symmetry can 
forbid lower order contribution to the neutrino mass. In sec. 4, we discuss how the Z n symmetry 
can be generalized to other groups. In sec. 5, we discuss the restrictions from LFV rare decays and 
anomalous magnetic moment of the muon. In sec. 6, we summarize our conclusions and briefly 
comment on implications for LHC signatures. 



2 General setting of the models 

We extend the Standard Model by introducing Ns scalars, Si, i = 1, . . . , Ns, and Np left-handed 
Weyl fermions, Fj, j = 1, . . . , Np. We assume that the leptons constitute a portal to the hidden 
sector via a Yukawa coupling of form 

N s Np / \ 

Cy = y~iy~l gijaSjFjLg, where L a = I "\ . (1) 
i j V a J 

Throughout the paper, we will adopt a two component notation and write all fields as left-handed 
Weyl spinors, i.e. in the (|, 0) representation. The product of two Weyl spinors £, rj is therefore 
defined by £77 = £ T cr/, where c denotes the charge conjugation matrix (i.e., en = C22 = 0; c\% = 
— C21 = 1). We will use an index- free notation, unless a special discussion of the Lorentz structure 
is required. Notice that only the combinations of form FjL a in Eq. (1) are allowed. A combination 
of form FjL a is forbidden by the Lorentz symmetry. We focus on the neutrino mass generation 
via coupling Eq. (1) and do not consider other radiative neutrino mass generation mechanisms, 
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e.g., via a coupling of the new particles to RH charged leptons (See e.g. [9]) or via two W boson 
exchange [15]. 

Of course, Cy should be a SU{2) x U(l) invariant combination. As a result, the sum of 
hypercharges of Si and Fj is opposite to the hypercharge of L a , i.e., Ys i + Yp j = —Yl = 1. 
In case that the hypercharge of the chiral fermions that we are adding is non-zero, anomaly 
cancellation might require addition of extra chiral fermions. There are various ways to make 
the combination invariant under SU{2). For example, if Fj is a triplet and Si is a doublet, the 
combination e a bSf (Fj) bc L c a is a SU{2) invariant (a, b and c are SU{2) indices). Our discussion 
of the implications of the Z n symmetry for loop contributions to the neutrino mass matrix is 
independent of the behavior of the fields under SU{2) so we shall not specify the behavior of the 
fields under the electroweak symmetry. 

Taking G v = Z n , the fields in Eq. (1) transform as follows 

S. e^^Si , 
Fj -> e^^Fj , 

L a — > L a , (2) 

such that 

a Si + &Fj 

as + «f = mod n 44> — G Z 44> as + «f £ nL. (3) 

n 

If we promote Z n to C/(l), the condition in Eq. (3) should be replaced by as^ + ap 4 = 0. 

Obviously, a Majorana mass term for the fermion F is forbidden, unless the Z n charge of F 
fulfils 2aF = mod n. A Weyl fermion F with 2a f ^ mod n therefore needs another Weyl 
fermion, F' with = — api to form a Dirac mass term. In case that F and the conjugate of F' are 
in the same representation of SU{2) x U(l), the mass term will be simply of form FF'. Anomaly 
cancellation in this framework will be automatic. If F and F' are in different representations of 
SU{2) x U(l), anomaly cancellation might require additional chiral fields and a mass term can 
emerge only after electroweak symmetry breaking. For example, if F is a doublet and F' is a 
singlet, the mass term can originate from a term of form FF'H. 



3 Loop contributions to neutrino masses 

As discussed earlier, we focus on models within which a Dirac neutrino mass is forbidden by a 
Z n symmetry and Majorana neutrino masses are produced only at loop level. In subsection 3.1, 
we make general remarks on the loop contributions to the neutrino mass. In subsection 3.2, we 
focus specifically on the one-loop contribution. In subsection 3.3, we discuss the conditions for 
constructing a lower loop contribution to the neutrino mass using the propagators and vertices in 
a general multi-loop diagram, with a specific discussion about the two-loop and three-loop cases 
in the subsequent subsections. We will analyze the different possible topologies without specifying 
the SM model charges or even the number of new fields. 
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3.1 General remarks on the loop level neutrino mass 

We consider diagrams contributing to neutrino masses and we will indicate the scalar propagators 
by dashed lines and the fermion propagators by solid ones. In general, the scalar propagator can 
involve more than one scalar: /iSi-S^y. Notice that a propagators of form (S1S2) can be rewritten 

in form of (^SiS^ by redefining S3 = S%- Without loss of generality, we will work in a basis with 

diagonal kinetic terms; as a result, a propagator of form (^FiF^ for F\ ^ F2 does not exist. The 
propagators of form (^FF^ preserve any Z n symmetry. In general, the fermionic propagator can 
be either chirality flipping (i.e., of form (FiF^c) or (cF^F^)) or chirality conserving (i.e., of form 
(FF*)). Propagators of type (fiff c) and (cF*F^ can result from Dirac or Majorana mass 
terms. 

In order to generate a Majorana mass term vv for neutrinos, lepton number has to be broken 
by two units. This can be achieved in various ways. An extensively studied option is to have 
a Majorana mass term for the new fermions and a lepton number violating mass term of form 
m 2 S 2 /2 for the new scalars. However, the options in general are wider. For example in case of the 
two-loop diagram in Fig. 5c, if we assign lepton number equal to —1 to S\ and S2 (or to Fi and 
F2), lepton number will be broken by two units by the S1S2S3 vertex or the -Fii^S^ vertex. 

A Majorana neutrino mass term vv can arise only after electroweak symmetry breaking. As we 
already briefly mentioned, the Weinberg operator is the lowest dimension operator that can induce 
Majorana mass for neutrinos. Higher dimension operators can also result in neutrino mass. See 
[4] for a classification of all AL = 2 operators leading to neutrino masses. However, additional SM 
fields, which are not charged under G v do not affect our discussion of the G v symmetry besides 
generating the neutrino mass at a higher loop order. A particular class of these operators are 
(H^ H) m (HL)(HL) with more than one pair of Higgs fields attached to the diagram, which have 
been studied in [5]. In the following, we will concentrate on the simplest origin and only discuss 
the Weinberg operator (HL)(HL), i.e., a pair of Higgs H being attached to the loop diagram 
giving mass to neutrinos. Let us discuss each option separately. 

• The Higgs can be attached to a vertex of type S1S2S3 via a renormalizable coupling HS1S2S3, 
provided that this combination forms a singlet. However, it cannot be attached to a fermionic 
vertex because the corresponding term in the Lagrangian will be non- renormalizable. 

• Let us now discuss the case in which the Higgs field is attached to the propagators. Propa- 
gators of form (FF*) and (S&) cannot break SU(2) x U(l) but propagators of type (SS) 
or (FF T c) in principle can do so. Let us take a general propagator of form (4>^) where cp 
and ip are either both scalars or both left-handed fermions. If ((f)ip^) breaks hypercharge by 
one (two) units, electric charge conservation implies T^((j)) — T^fy) = 1 (2). This means that 
4> and ip cannot be both singlets. A propagator involving only one Dirac field in an SU(2) 
doublet representation can break hypercharge by at most one unit. However, if we allow 
more than one field to be involved in the propagators, more possibilities open up. The line 
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Figure 1: Fermion propagator breaking the hyper-charge by two units 

shown in Fig. 1 is an example. To have such a line, the required terms are HF3F2, HF^F\ 
and mF 3 F^F^. From now on, for brevity we shall not emphasize on the requirement of mass 
insertion for a line such as the one denoted by F3 and F3 in Fig. 1. 

Consider a loop that contains nj internal lines plus nsv vertices that involve three scalars. For 
such a diagram, there are 

ni [(ni + l)/2 + n S v] + n sv (n S v ~ l)/2 

ways to attach the pair of external Higgs fields to the diagram. For example, in case of Fig. 6b 
ni = 8 and nsv = 3 which means there are 63 ways of attaching the external pair of Higgs fields. 
To avoid cluttering the figures with this plethora of possibilities, we do not show the Higgses 
attached to the diagrams. 

3.2 One loop 

At one loop, there are two possible diagrams, which are not accompanied by a tree-level contribu- 
tion. They are shown in Fig. 2. Had we included neutral fermions invariant under G v or allowed 
new scalars to develop a VEV, we could have more types of one-loop diagrams [6]. 

The propagators in the one-loop diagram Fig. 2a are of form {S1S2) and {F\F^c) where S2 
and i*2 may or may not be the same as £1 and F\, respectively. Let us denote the Z n charge of an 
arbitrary field <f> by a^. In order for the propagators to be Z n invariant, the charges of Si and S2 
as well as the ones of F\ and F2 have to add up to mod n. The existence of the vertices leads 
to similar conditions. The following relations need to be satisfied 

asi + as 2 > a F! + ccf 2 , &Si + «5 2 + a F 2 ^ n ^ 1 (4) 

which lead to 

as ± = — as 2 = a F 2 — ~ a F 1 mod n . (5) 

In the specific case that the pair of (Si,F\) is identified with (£2,-^2), we find 2as 1 = 2ai?i — 
ols\ + a Fi =0 mod n and if there is no other field 4> with non-trivial Z n parity, any choice for n 
will be equivalent to n = 2. However, in general when S\ 7^ S2 and F\ 7^ F2, n might be different 
from 2. 
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(a) One loop diagram with mass 
term SS on scalar line 



(b) One loop diagram without 
mass term SS on scalar line. 



Figure 2: Effective neutrino mass generation at one loop. 




Figure 3: General examples of loops contributing to the neutrino mass which can be disconnected 
by cutting a pair of propagators. 

Let us discuss the special case S = Si = S%, which corresponds to the diagram shown in 
Fig. 2b. Apparently, in this case the scalar line, (SS^) cannot break SU(2) x U(l) which means 
both Higgs fields have to be attached to the fermion line. As discussed earlier, if the fermion line 
is needed to break hypercharge by two units, it has to involve at least one fermion in addition 
to F\ and F% (see Fig. 1). Notice that in this case no lepton number violating mass insertion of 
type m 2 S 2 /2 is required. Instead, the simultaneous presence of the SF\L and S^F^L vertices and 
the F\F2 mass term breaks lepton number. The Lorentz structure of the Majorana mass term, 
vis, cannot be created by a fermion propagator of type (FF^) so a mass insertion of the fermionic 
propagator is required. 

The simplest example of this type of models is Ma's scotogenic model [11], which introduces 
one additional inert Higgs doublet rj as well as three RH neutrinos N{. Neutrino masses are 
generated with the internal propagators (NN T c) and (ry°r/°) with rf being the neutral component 
of the inert Higgs doublet rj. 

The one loop suppression does not suffice to explain the smallness of neutrino masses by 
itself. A further suppression is needed, which might be due to the smallness of the lepton number 
violating mass insertions that are indicated by crosses compared to the overall masses of the 
particles propagating in the loops. Another explanation of the additional suppression might be a 
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sequence of symmetry breaking which naturally suppresses certain couplings (See e.g. [14]). 
3.3 Reduction of multi-loop contribution to one-loop 

Let us consider a general loop contributing to the neutrino mass which can be disconnected by 
cutting a pair of propagators as shown in Fig. 3a. If, as shown in Fig. 3a, the fermionic propagator 
is of chirality flipping nature (FiF^c), the vertices of the following types will be allowed by the 
G u = Z n or £7(1) symmetry: 

S\L a Fi and S 2 LpF2 . (6) 

If both the Higgs fields are attached to this pair of lines, these two vertices can be made S£7 (2) x 
£7(1 ) invariant, too. As a result, a one- loop diagram contributing to the neutrino mass with the 
couplings in Eq. (6) can be formed. 

On the other hand, if the propagator is of the chirality-flipping form (cF^F^), the vertices 
S\L a Fi and S\l^F 2 are allowed by the G v = Z n or £7(1) symmetry but they may violate £7(1) 
hypercharge. Thus, unlike the previous case, the presence of a one-loop contribution to the neutrino 
mass is not guaranteed. On the other hand, if the fermionic line is of the chirality-conserving form, 
i.e., (FF^), see Fig. 3b, vertices in Eq. (6) might violate the G u = Z n or £7(1) symmetry and 
again, a one-loop contribution to the neutrino mass does not necessarily exist. 

In summary, if there exists a multi-loop contribution to the Weinberg operator, (HL a )(HLp) 
compatible with G u = Z n or £7(1) which can be disconnected by cutting a pair of fermionic and 
scalar lines, there will be also a one-loop contribution to the neutrino mass provided that (i) both 
Higgs fields are attached to these two lines; (ii) the fermionic propagator in question is chirality- 
flipping and of form (FiF^c). 

Let us now consider a general multi-loop diagram of form shown in Fig. 4a in which the internal 
loop only gives a correction to the wave function of the scalar. Topologically such a diagram is 
distinguished from the rest by the fact that by cutting the scalar lines directly connected to u a 
and up, this line will be disconnected. Figs 5a, 5b, and 5c are examples of such diagrams but Figs. 
5e and 5f are not. If a contribution of this type exists, the G v = Z n or £7(1) symmetry allows 
a term such as S\S%, too. Depending on the electroweak behavior of S\ and S 2 , this mass term 
can result from terms such as e a ),SfH b S2 (for singlet S2 and doublet Si with Ys t + Ys 2 = — 1), 
e ah SlS b 2 (for doublets Si and S 2 with Y Sl + Y S2 = 0), e ab e cd SfH b S%H d (for doublets Si and S 2 
with Ys 1 + Ys 2 = —2) or S1S2 (for singlets Si and S2 with Y$ 1 + Y$ 2 = 0). As a result, these 
diagrams are always accompanied by a one-loop diagram. 

Let us now discuss diagrams of type in Fig. 4b in which the internal loop gives correction to the 
wave function of the propagating fermion. Similarly to the correction to the wave function of the 
scalar propagator, the wave function correction to the fermion propagator can also be written as 
a F1F2 mass term which respects the G u = Z n or £7(1) symmetry. However, if F\ and F 2 are both 
electroweak doublets, F1F2 will form an electroweak triplet. Thus, two factors of (H) are needed 
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(a) scalar line (b) fcrmion line 

Figure 4: Wavefunction renormalization of internal propagators. 

to contract it to a SU(2) x U(l) singlet. In other words, the corresponding term will be non- 
renormalizable. Hence, this diagram is not necessarily accompanied by a lower loop contribution 
depending on the electroweak structure of the fermions. 

Let us suppose a coupling of form g a SFL a compatible with G u = Z n or U(l) exists. There 
must be another F' with cx.pi = —ap to obtain a Dirac mass term for F (either directly or after 
electroweak symmetry breaking). The G v = Z n or U(l) symmetry does not forbid a term of form 
g' a S' F' L a . The neutrinos then obtain a Majorana mass proportional to g a g'n at one loop. In 
the discussion of higher loops, we implicitly assume that such a possibility is forbidden by other 
symmetries such as the electroweak symmetry. 

3.4 Two loop 

At two-loop level, there are more possible diagrams. In the following, we again take G u = 
Z n or U(\) and discuss in which cases the symmetries forbid the lower loop contribution to the 
neutrino mass. In Fig. 5, for the sake of simplicity, some of the scalar lines are marked by a single 
letter such as S\ and S^- If the external Higgs is attached to any of these lines, this line will in 
fact involve more than a single field. When the fermionic propagator involves an even number of 
fields, the chirality will be flipped because it requires a mass term such as F%F2- In the diagrams, 
the arrow indicates the direction of the flow of the G u charge, both for fermion propagator and 
scalar propagator. Notice that for the fermion lines, the chirality might flip but the direction of 
the arrow will remain the same. If a chirality flip is required (e.g., Fig. 5f), the fermionic line 
will involve the fermion, F and its partner F' with opposite G u charge that together form a Dirac 
mass. This simplified way of marking does not generally affect our discussion below. We will be 
more specific when it does. 

As discussed in the previous section, the diagrams 5a, 5b, and 5c are always accompanied 
by a one-loop diagram but this is not necessarily the case for the diagram in Fig. 5d, as it has 
been discussed in the previous section. The diagrams 5e and 5f cannot be reduced to a one loop 
diagram, as long as S3 and F3 transform non-trivially under G v = Z n or U(l), so they can give 
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Figure 5: Two loop diagrams 

the dominant contribution to neutrino mass. 

The two-loop suppression factor (g 2 /(167t 2 )) 2 Mjvevk for m NEW ~ 100 GeV and g ~ 0.1 will 
not be enough to explain the smallness of the neutrino mass but additional sources of suppression 
are necessary. For example, like the case of the Zee-Babu model [3], the suppression can result 
from the heaviness of S3 (P3) in Fig. 5e (Fig. 5f). In the following, we are going to discuss the 
possibility of a further loop suppression, when neutrino masses are generated at three loop order. 

3.5 Three loop 

Before starting the discussion of the three-loop diagrams, we emphasize that the comment in the 
first paragraph of the previous section on marking the propagators applies here, too. The three- 
loop diagrams contributing to the neutrino mass can be divided into three categories: (i) diagrams 
in which the inner loops correct the wave function of the internal lines. Such diagrams are already 
discussed in sect. 3.3. (ii) Planar diagrams shown in Figs. 6, 7, 8 and 9. We shall discuss all 
these diagrams in detail below in context of G u = Z n and G v = U(l). (iii) Non-planar diagrams 
shown in Fig. 10 which will be discussed in detail in this section in context of G u = Z n and briefly 
for G v = U(l). Although at first sight it seems there are four non-planar three- loop diagrams 
contributing to the neutrino mass, only two are independent. As demonstrated in Fig. 10 twisting 
the vertices denoted by Pi and P3, the diagrams on the left- and right-hand sides of Fig. 10 convert 
into each other. There are therefore only two distinct non-planar diagrams. 
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(a) 



(1)) 



(c) 



5*i 



*3\ 



i 5*5 45 2 




(d) 



Figure 6: Planar three loop diagrams. 



The G v = Z n or G v = U(l) symmetries do not forbid lower order loop contributions for any of 
the planar diagrams besides the one in Fig. 9, but the pattern of electroweak symmetry breaking 
as well as the requirement of chirality flipping might prevent some. Let us discuss this possibility 
in detail. In Fig. 6 and Fig. 7, the internal loop in the red dashed dotted circle can be replaced by 
an effective G v conserving vertex. Let us first consider the diagrams in Fig. 6. For a Lagrangian 
symmetric under G v = Z n (or under G v = U(l)), the presence of these diagrams is possible only 
if aSi associated with Si, S2 and S3 add up to an integer times n (or add up to zero). Thus, 
the G v invariant Lagrangian contains vertices of type S1S2S3 + h.c. unless it is forbidden by some 
other symmetry. In particular let's consider the SU(2) x U(l) symmetry: each line might involve 
more than one field with different SU(2) x £7(1) quantum numbers and we should specify the 
fields that directly leave the red dashed dotted circle. In case that Y$ t + Y$ 2 + Y$ 3 = or —1, 
the corresponding vertex can be just S1S2S3 + h.c. or 51525*3 (H) + h.c, respectively. However, 
there is no renormalizable vertex of form S1S2S3 if Ys 1 + Y$ 2 + Y$ 3 < — 1 and no corresponding 
two-loop diagram. Similarly, in case of the diagrams in Fig. 7, the loop in the dashed dotted 
red circle can be also replaced by a renormalizable Z n invariant vertex of form S\f\F2 unless 
Ys 1 + Yp 1 + Yp 2 7^ 0. Notice that we implicitly assume that both left-handed fields F\ and F2 
leave the red dashed dotted circle; i.e., the corresponding effective vertex is of form S^F^F^ rather 
than Sj-Frj F\ or S\f'[F2 which are forbidden by Lorentz structure. Let us consider the loop on the 
right-hand side in the diagrams of Fig. 8. If, as indicated in the figures, the fermion marked with 
F enters this loop, the loop can be replaced by a renormalizable Z n invariant Yukawa coupling of 
form S*FL/3, which conserves hypercharge if Ys = Yp — 1. 
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Figure 7: Planar three loop diagrams. 




(a) (b) (c) 

Figure 8: Planar three loop diagrams. 

The implication of Z n for Fig. 9 as well as the non-planar diagrams in Fig. 10 is more compli- 
cated. In particular, Z n does not always allow them to be accompanied by a dominating two loop 
contribution to the neutrino mass. Let us first consider diagram 9. The Z n symmetry implies 

a Fx = —&Sx + a F 3 = 0(s 3 — CtSx + n &3 (7) 

a>F 2 = —as 2 + nk 2 as 4 = asx + «s 2 — a s 3 + 

where k{ are arbitrary integers. We are interested to find out whether there is a Z n symmetry with 
certain assignment which is compatible with the diagram in Fig. 9, i.e., satisfies Eqs. (7) but 
forbids lower loop contributions. To answer this question, we have solved equations (7) under the 
condition that none of the one- and two-loop diagrams respectively in Fig. 2 and Fig. 5, is allowed 
by the Z n symmetry. The values of and opv can be set such that the Z n symmetry forbids 
the lower orders of contributions to the neutrino mass for n > 16. One example is the charge 
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Figure 9: Planar three loop diagrams. 




assignment 

= 1 «5 2 = 3 as 3 = 9 as 4 = 11 ap 1 = 15 ap 2 = 13 ap 3 = 8 (8) 

for Ziq. However, for smaller values of n such an assignment is not possible. 

Similarly in case of the non-planar diagrams, for certain assignments of ag. and op., the Z n 
symmetry forbids all diagrams with a lower loop order. Let us consider diagram Fig. 10a. The Z n 
symmetry leads to the following relations for the different fields propagating inside the loops: 

ap 1 = —ds 1 + nk\ as 3 = as 5 — as 2 + nk^ (9) 

aj? 2 = — as 2 + n&)2 Q-Sa = a Si + a s 2 ~ a s 5 + 

ap 3 = as 1 — as 5 + nk 3 . 

We have found that the smallest value of n for which the Z n symmetry forbids lower order loop 
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Fig. 10a Si S 2 S 3 S 4 S 5 Fi F 2 F 3 



Fig. 10b F 2 F 1 F 5 F A F 3 S 2 S t S 3 



Table 1: Translation table between Fig. 10a and Fig. 10b. 



contribution to the neutrino mass is n = 16. One particular example is 

«5 1 = 2 as 2 = 6 as 3 = 13 as 4 = 5 as 5 =3 ap 1 = 14 ap 2 = 10 ap 3 = 15 (10) 

for Ziq. That is for n < 16 any possible assignment of and will also lead to a dominant 
lower loop diagram. In this analysis, only the Z n symmetry is considered. Of course, depending 
on the field content of the specific model, the pattern of the hypercharge breaking as well as the 
form of chirality flipping might also forbid the lower loop contribution. 

A similar analysis is performed for the diagram in Fig. 10b. To be compatible with the Z n 
symmetry, the following set of equations has to be satisfied 



for any integers k{. Here, n = 16 is again the smallest value of n for which Z n forbids a lower loop 
contribution to the neutrino mass. One example of charge assignment for Z%% forbidding lower 
loop contributions is 

as 1 = 2 as 2 = 6 as 3 =3 ap x = 14 ap 2 = 10 a^ 3 = 7 «f 4 = 1 a^ 5 = 9 . (12) 

Notice that all three of the three loop diagrams require n ^ 16 to forbid lower order loop diagrams. 
We will first show that the conditions for Fig. 10a and Fig. 10b are equivalent. If we identify the 
charges in Eq. (9), as indicated in Tab. 1, it is straightforward to show, that we obtain a set of 
equations, which is equivalent to Eq. (11) and vice versa. Note that we are identifying the charges 
of the fermions in one diagram with the ones of the scalars in the other diagram and vice versa. 
Remember that to form a Dirac mass term for the new fermions, we have to introduce a partner 
for each fermion so the charges of scalars and fermions may be treated on an equal footing. As a 
result, the equivalence of the conditions for the absence of lower loop contributions directly follows. 
The equivalence of charge assignments also implies that for a given n ^ 16, there must be exactly 
the same number of possible charge configurations for each diagram that forbids the lower order 
loops. In fact, solving the equations for n = 16, we find that there are exactly 32 possible charge 
configurations for each of the diagrams Fig. 10a and Fig. 10b that forbid lower order contributions. 

It is straightforward to verify that after replacing 



a.p 1 = —ots 1 + nki 
ap 2 = — as 2 + nk 2 
aF 3 = —cts 2 ~ a s 3 + n ^3 



ap 4 = — as 1 + as 3 + nk^ 
«F 5 = a Sl - a S2 - a S3 + nk 5 



(11) 



as 3 o a Ss 



CtF 3 — CXp z 



3 ) 



(13) 
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Eqs. (7) will be equivalent to the first four equations in Eq. (9). As a result, if for a given n there is 
an assignment of charges for the field content of Fig. 10a that forbids lower loop contributions to the 
neutrino mass, the corresponding assignment for Fig. 9 will also forbid lower order contributions. 
However, the opposite statement is not valid. In fact, the presence of an extra scalar field in case 
of Fig. 10a gives more freedom to construct more possible lower loop diagrams. Setting n = 16, 
we find 112 solutions for Fig. 9 that forbid lower loop contributions but only 32 such solutions for 
diagrams in Fig. 10a. 

In this discussion, we have assumed that each fermionic (scalar) propagator has an independent 
a Fi ( a Sj) value. The relation between and comes from the requirement that this diagram 
respects Z n . However, within specific models, there might be more restrictions. For example, 
let us assume that S\ and 52 in diagram Fig. 10a are the same fields. We then conclude that 
for any value of n and any choice of ap i and 0:5. (provided that ots 1 = otg 2 ), the Z n symmetry 
allows the diagram in Fig. 10a to be accompanied with two diagrams of form in Fig. 5c with 
the following replacements of the fields in Fig. 5e: (Si, 5*2, S3, Fx, F%) — > (Si, S3, S5, Fx, F%) or 
(Si, S2, S3, Fx, F2) (Si, S3, S4, Fx, -F3) where F3 is the Weyl fermion with charge opposite to 
that of F3 that together form a Dirac mass term F3F3. Depending on the field content of the 
model, it is possible that one or both of these two-loop diagrams are forbidden by SU(2) x U(l) or 
by pattern of chirality flipping. For example, if in Fig. 10a, one of the Higgs fields is attached to S4 
and the other to S5 , neither of these two- loop diagrams can exist because each violates hypercharge 
only by one unit. 

Let us now briefly discuss the possibility of replacing Z n with a continuous U(l). In this 
case setting n = 0, relations in Eqs. (9) and (11) remain valid. For a U(l) with general ag ; and 
OLFi the three- loop non-planar diagrams and the planar diagram in Fig. 9 cannot be accompanied 
with a lower loop contribution unless in very specific cases. For example, in the specific case that 
as 4 = — as 5 , there might be also a one-loop contribution to the neutrino mass accompanying the 
diagram in Fig. 10a. 

4 General symmetry 

In previous sections, we focused on the implications of G v = Z n or symmetry on the neutrino 
mass generation at loop levels. As we shall discuss below, some of these discussions can be applied 
for a general symmetry group G v . We will make a further generalization in this section. Motivated 
by the DM models, in the previous sections, we have assumed that the SM particles are all invariant 
under the G v symmetry. In this section, we discuss the consequences of relaxing this assumption. 
We will still assume that none of the scalar fields which non-trivially transform under G u (i.e., are 
not invariant under G v ) receives a VEV. Thus, the SM Higgs is invariant under G v . 

Let us reconsider coupling (1) and review the results that we found in the previous sections. 
• As long as all new neutral fermions non-trivially transform under G v , the Dirac mass term 
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for neutrinos will be forbidden by this symmetry to all orders in perturbation theory. 

Consider a general n-loop diagram contributing to neutrino mass. Suppose that there is a 
sub-diagram within this diagram that absorbs scalar lines Si, S2 and S3. The existence of 
such a sub-diagram implies that a combination of form S1S2S3 is invariant under G v . If 
no Higgs VEV is attached to this sub-diagram, it means this combination is also invariant 
under electroweak symmetry and this term is allowed in the Lagrangian. Thus, there should 
be a lower order contribution to the neutrino mass where this sub-diagram is replaced by 
the S1S2S3 vertex. Similarly, if there is only one Higgs VEV attached to this sub-diagram, a 
renormalizable coupling of form S1S2S3H exists in the Lagrangian which can lead to a lower 
order contribution to the neutrino mass. However, if both Higgs VEVs of the Weinberg 
operator are attached to the sub-diagram, the corresponding electroweak vertex will not be 
renormalizable. Similarly, there might be a sub-diagram in which two fermions (Fi and 
F2) and one scalar (Si) leave the sub-diagram. Examples are shown in Fig. 7, where the 
sub-diagram is inside the red circle. Regardless of the details of G u and whether the SM 
particles are invariant under it or not, a vertex of form F1F2S1 is invariant under G v . If 
it is also invariant under electroweak symmetry (i.e., if no Higgs VEV is attached to the 
sub-diagram) this renormalizable term will be present in the Lagrangian. This means this 
diagram is accompanied by a lower order diagram in which the sub-diagram is replaced by 
vertex of form F1F2S1. 

Let us now consider a general diagram that contains a sub-diagram which is a correction to 
the self-energy of a scalar line, that is the external lines that are attached to this sub-diagram 
are two scalars Si and S2 (See Fig. 4a). The S1S2 term will be G u invariant and in case that 
less than three Higgs VEVs are attached to this sub-diagram, it can be made electroweak 
invariant by adding appropriate number of Higgs fields (i.e., S1S2, S1S2H or S1S2HH). 
The diagram is accompanied by another one in which the sub- diagram is replaced by the 
corresponding renormalizable vertex. Similar consideration holds for a sub-diagram giving 
correction to the fermion self energy but in this case more than one Higgs fields cannot 
be added otherwise the corresponding term will be non-renormalizable. Notice that this 
consideration holds valid regardless of the details of the G u symmetry and the behavior of 
the SM fermions under this symmetry. 

The theorem of section 3.3, regarding diagrams of form shown in Fig. 3 holds valid for a gen- 
eral G v independent of the behavior of the SM fermion under the symmetry transformation. 

As shown for the special case of G v = Z n , the G v symmetry can forbid lower order con- 
tribution for the two-loop diagram of topology shown in Fig. 5e and Fig. 5f, as well as the 
three-loop diagrams of topology in Figs. 9 and 10. 
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Figure 11: LFV diagram. The photon is emitted either from the initial or final state or from the 
charged particle in the loop. 

5 Lepton flavor violation 



After discussing the general structure of the different diagrams leading to neutrino masses, we 
will discuss lepton flavor violating rare decays which have proven to lead to strong constraints on 
radiative neutrino mass models. The general Yukawa coupling in Eq. (1) includes LFV couplings 
of left-handed charged leptons, l~, gij a Sifjl~ where Sj and fj are components of the multiplets Si 
and Fj such that the following relation holds valid between their electric charges: 

Qs i + Q fj = -Q l - = i. 

This LFV coupling will lead to LFV rare decay l~ — > + 7 as shown in Fig. 1 1 . Note that in 
contrast to the case of contributions to the neutrino masses, the LFV rare decays are generically 
allowed at one loop order unless there is a flavor symmetry forbidding the one loop contribution; 
see e.g. [16]. Neglecting the corrections of order of (m/„/m; a ) 2 , the decay rate is [17] 
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Notice that this relation holds regardless of the representation of the electroweak symmetry to 
which the new particles belong. Larger representations can lead to several different possible loop 
diagrams. 



16 



In the case that the neutrino mass originates from the three-loop contribution, the coupling 
should be of order one for m^EW ~ 100 GeV to account for m u ~ \J Am^ tm . This will lead to 
Br(/i — > erf) exceeding the present bound. However, there are ways to avoid these bounds making 
use of a particular flavor structure. In the following, we explain a simple and natural solution. To 
reproduce two nonzero neutrino mass eigenvalues, more than one pair of (S, F) coupled to L is 
required, which we will call (Si, Fx) and (£2, -^2)- Let us suppose that Fx only couples to L e while 
F2 couples to and L T . That is gxx/j, = 9xir = 922e = 0. In this case, the couplings conserve L e 
so — > e7 will be absent. As it is well know [18], the conservation of L e leads to a vanishing first 
row and column of the neutrino mass matrix, i.e., m ea = for a = e, fj,, r, and can therefore only 
serve as a leading order approximation of neutrino masses. If L e is softly broken by trilinear scalar 
couplings, the vanishing elements can be reproduced and the observed neutrino mass pattern can 
be reconstructed. Obviously, the breaking introduces fi — > &y, but it is controlled by the smallness 
of the symmetry breaking and will not be dangerous. 

The new coupling can lead to a new contribution to anomalous magnetic moment of muon as 
follows 

(g ~ g)g _ ml ^ | gtJ> | 2 

2 192vr 2 ^ ml 

ij s * 

Notice that for ms,rriF ~ 100 GeV and g ~ 0.5, this contribution can explain the observed 
anomaly. 



mi 



mi 
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mi 



mi 



6 Conclusions and discussion 

The smallness of neutrino masses is one of the longstanding problems in the phenomenology of 
particle physics. Various models have been proposed in the literature within which neutrino masses 
are produced at loop level so its smallness is natural. In this paper, we have discussed the class of 
models within which neutrino masses are produced at loop level via a Yukawa term that couples 
neutrinos to new scalars and fermions. We studied and outlined some general results that can be 
drawn from the topology of neutrino mass diagrams or the symmetry, G v , imposed on the model. 
We have discussed conditions on the G u symmetry and topology of loop diagrams that forbid the 
presence of a lower order, and consequently dominating, contribution to the neutrino mass. Under 
these conditions these diagrams will therefore give the dominant contribution to the neutrino mass. 
More general results are outlined item by item in sect. 4. 

In this paper, we have assumed that the G v symmetry remains unbroken. In case that the 
SM particles are invariant under the G u symmetry, the lightest new particle with a non-trivial 
behavior under this symmetry will be stable. If this particle is neutral, it can contribute to the 
dark matter in the universe. Independently of a given model, there have been studies of the impact 
of different symmetries on the DM predictions (See e.g. [8]). We shall discuss this issue in more 
detail elsewhere. 
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Within the class of models that we have discussed in this paper, the scale of new physics can 
be as low as the electroweak scale. The new particles that are added have no strong interactions; 
however, they can have electroweak interactions. The model can include new charged particles 
coupled to the Higgs which, along with the SM contributions, may explain the possible excess in 
the diphoton Higgs decay channel [19]. 

At a hadron collider such as the LHC, the only production mode of these particles is through 
electroweak interactions but in a lepton collider such as the ILC, these particles can be also 
produced via Yukawa interactions in Eq. (1) in the i-channel. If all the SM particles transform 
trivially under G u and all new particles carry G u charges, these particles can be produced at 
colliders only in association with other new particles such that the final products form a singlet 
of G u . For example, if G v = Z2, the new particles that are odd under Z2 can be produced only 
in pairs. Moreover, the G u symmetry implies that the decay products of new particles include 
lighter new (beyond SM) particles. In fact these particles will go through a chain of successive 
decays until they produce stable new particle(s). If these final stable products are neutral, they 
will appear as missing energy but, if they are charged, they can be detected. Since they are heavy 
and have no strong interactions, they will generally lose energy with a rate smaller than the muon 
energy loss rate which means they come to rest only after they exit the detector. In this case, the 
signature of the model will be quite distinct from the SM background, raising the discovery chance 
of the model. 

Decay of the new particles can take place through the Yukawa coupling in Eq. (1) which 
means each decay produces a lepton along with a new particle. The branching ratio to different 
flavors is determined by gij a . This is the same coupling that determines the flavor structure of 
the neutrino mass matrix. In principle, by studying the flavor composition of the decay products 
of the new particles, one can cross check these models. This possibility has been studied in detail 
in [20] for the specific case of the SLIM model [12, 13] where G v = Z2 and the new fermions are 
neutral. Unfortunately, the high rate of background and the uncertainty in luminosity will limit the 
capability of the LHC to extract the flavor structure of the coupling. However, a lepton collider 
can have a better chance of determining the coupling. If the model contains multiply charged 
particles, their production will be enhanced by square of their charge. Moreover their successive 
decay to multiple charged leptons plus missing energy or new stable charged particle will provide 
a distinct signature enhancing the chances of discovery at lepton and hadron colliders. 

One of the key ingredients used in our setup to suppress lower loop contributions is the 
presence of discrete symmetries at quantum level. Although in our setup discrete symmetries 
are global and not gauged, there are arguments that all symmetries, including discrete global 
symmetries, should be gauged in a theory of quantum gravity [ ]. This is due to the expectation, 
mainly advocated in string theory, that all symmetries have a geometric origin and the space- 
time itself is a locally constructed, secondary concept, in these settings. Thus, all the symmetries 
should be local symmetries. These symmetries could be broken due to quantum gravity effects 
and/or anomalies [22]. The effects of symmetry breaking (if there is any) are expected to be 
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suppressed by inverse powers of mp\. Based on the dimensional analysis, we can estimate the 
contribution from the G u violating effects to the neutrino mass to be mNEwi^NEW /fnpi) 2 ~ 
1CT 6 eV(mNEW /100 GeV) n+1 . Thus, we can safely neglect this effect. 

In summary, small neutrino masses can be generated at the loop level in models in which the 
leptons couple to the new sector. New symmetries guarantee the absence of a Dirac mass term 
for neutrinos and can forbid lower loop diagrams. The presence of this symmetry leads to neutral 
stable candidates which might explain the observed dark matter abundance in the Universe. The 
additional suppression due to the higher loop order allows to lower the scale of new physics down 
to the TeV scale keeping large couplings and providing specific testable signatures at colliders and 
observable lepton flavor violating processes. 
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